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Effects of pressure gradients 
on turbulent premixed flames 

By D. Veynante 1 AND T. Poinsot 2 


1. Motivations and objectives 

Most ducted turbulent flames are subjected to external pressure gradients. Com- 
pared to “free” flames, i.e. turbulent flames without externally imposed pressure 
gradients, the combination of the external pressure gradients with the large density 
changes found in premixed flames may lead to strong modifications of the flame 
structure. These modifications are mainly due to the differential buoyancy effects 
between cold, heavy reactants and hot, light products. They affect turbulent trans- 
port along with many characteristics of the flame itself, such as the flame speed, 
thickness, wrinkling, and local structure. Pressure gradients are also a key mech- 
anism for the counter-gradient turbulent transport described below. Accordingly, 
studying the effects of pressure gradients on premixed turbulent flames is an im- 
portant issue both for fundamental understanding of turbulent combustion and for 
modeling. 

Using the assumption of single-step chemistry, the mass fractions of the reactive 
species are all linearly related (Williams 1985) and may be expressed in terms of a 
single reduced mass fraction: the reaction progress variable c. The progress variable 
ranges from zero to unity in the fresh and fully burnt gases, respectively. Using the 
classical Favre decomposition, a quantity q can be split into a mass-weighted mean, 
q = pq/p, and a turbulent fluctuation, q" . The transport equation for the mean 
reaction progress variable c may be written as: 

dpc dpu x Z dpu”c n _ dJk — 

dt + dxi dxi dxk 

where p is the mass density, U{ is the flow velocity, Jk is the molecular diffusion 
flux, Cj c is the volumetric production rate of the chemical reaction, and the over- 
bar superscript denotes conventional Reynolds ensemble-averaging. Equation (1) 
has the form of a standard turbulent transport equation where the rate of change 
of c results from a balance between convection by the mean flow, convection by 
the turbulent flow, molecular diffusion, and chemical reaction. The contribution of 
molecular diffusion is usually neglected for high Reynolds number flows. In Eq. (1), 
two terms need to be modeled: the mean reaction rate u> c and the turbulent trans- 
port pu"c n terms. The first term has received considerable attention in recent years 
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and various models have been derived and incorporated into practical codes for 
turbulent combustion. The second term, however, has received considerably less 
attention and is generally described with a simple classical gradient eddy-viscosity 
model: 

Ht oc 


pu'fc" =pu'fc" = - — — 


( 2 ) 


where nt denotes the turbulent dynamic viscosity and <7 C a turbulent Schmidt num- 
ber. 

Both theoretical and experimental research (Bray et al., 1981, 1982; Shepherd et 
al., 1982) have shown the occurrence of c ounter -gradient transport in some turbulent 
flames: flames where the turbulent flux pu”c" and the c gradient, dc/dxi, have the 
same sign in opposition with the prediction of Eq. (2). This is generally due to the 
differential effect of pressure gradients on cold reactants and hot products. Recent 
studies based on direct numerical simulations of turbulent premixed flames without 
externally imposed pressure gradients (Trouve et al., 1994; Rutland & Cant, 1994) 
have confirmed that counter-gradient diffusion was found in simulations, but that 
classical gradient diffusion was also possible. A criterion indicating the presence of 
gradient or counter-gradient diffusion in atmospheric flames has been derived by 
Veynante et al. (1995). This criterion leads to a reduced number called Nb, or the 
Bray number: 


N b = 


2a(u' /s ° t ) 


( 3 ) 


where s° is the laminar flame speed, it' is the RMS turbulent velocity, r is the 
heat release factor defined as r = T h /T u - 1, with T being the temperature and 
indices u and b referring to the fresh and burnt gases, respectively. The term a is 
an efficiency function of order unity, introduced to take into account the reduced 
ability of small turbulent vortices to affect the flame front. For low values of Nb, 
typically Nb < 1 , where flames exist in relatively large turbulence intensity, gradient 
diffusion is obtained. For large values of Nb counter-gradient diffusion occurs. In 
fact, counter-gradient turbulent diffusion is promoted by heat release and thermal 
expansion (increasing values of r) whereas increased turbulence intensity tends to 
induce gradient transport. 

The work of Veynante et a/.(1995) was performed for free flames without exter- 
nally imposed pressure gradients or volume forces, such as gravity. Since turbulent 
transport in flames appears to be controlled by a dynamic balance between fresh and 
burnt gases, confined flames subjected to strong pressure gradients should exhibit 
a large sensitivity to these gradients. For example, we expect that imposing a pres- 
sure gradient on a turbulent flame exhibiting counter-gradient diffusion may lead to 
a gradient-diffusion situation. This change could affect the flame brush thickness, 
the turbulent flame speed, and finally the complete structure of the turbulent flame 
brush, as described in a number of papers by Masuya & Libby, 1981; Bray et al., 
1982. 

Our objective in this study is to explore the effects of pressure gradients on 
premixed turbulent flames using direct numerical simulations. We will first recall 
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the physics of this phenomenon and how pressure gradients may be included in a 
simulation for reacting flows in Section 2. We will also describe the most general 
theory to treat the problem, i.e. the Bray-Moss-Libby formulation. The simulation 
used for this work will be described in Section 3. Section 4 will present the structure 
of laminar flames submitted to pressure gradients. Section 5 will present results for 
turbulent flames. Finally, Section 6 will describe a model for the turbulent flux 
incorporating the effects of pressure gradients. 

2. Pressure gradients in premixed flames 

2.1 The Bray-Moss-Libby approach for turbulent transport in premixed flames 

Bray et al. (1989) have proposed a simple algebraic closure for the reaction term 
u>c in Eq. (1), but focus their attention on the turbulent transport term pu"c". In 
the Bray-Moss-Libby (BML) approach, the flame is analyzed as a thin flame sheet, 
or “flamelet,” separating fresh reactants (c = 0) and fully burnt products (c = 1). 
This assumption leads to a bimodal probability density for the progress variable c, 
and the turbulent flux is then expressed, according to Bray 1980, as: 


pu"c" = pc( 1 - c)(ui b - u iu ) (4) 

where W t% and xT tb are the conditional mean velocities within the unburnt and burnt 
gases, respectively. The occurrence of counter-gradient transport may be easily 
explained from this expression. Let us consider a left-traveling flame along the re- 
direction ( dc/dxi > 0). Thermal expansion and the associated flow acceleration 
through the flame, along with favorable buoyancy and/or pressure gradients, will 
tend to make ui b larger than i77 u , thereby promoting counter-gradient turbulent 
diffusion of c, resulting in pu"c" > 0. Under the Bray-Moss-Libby approach, the 
second and third turbulent moments such as pu"u ", pu"u"c", and ~pu”u"u'j c may be 
directly expressed as functions of conditioned quantities in fresh and in burnt gases. 
Nevertheless, conditional quantities such as ( ui b — tZ7 u ) are difficult to close and 
an alternative approach must be pursued for estimating the turbulent transport: 
pu"c". A simple algebraic closure based on the eddy viscosity concept cannot 
be used here . In t he Bray-Moss-Libby model, closure is achieved by a transport 
equation for pu'/c" (Bray 1980, Bray et al. 1989, Bray 1990). A brief derivation of 
this equation is provided here, including a constant volume force Ff and a constant 
acceleration T , . We start from the momentum equation: 


dpui 

dt 


dpuiUj 

dxj 


3P_ 

dii 


+ Ff + pTi + 


drjk 

dxk 


( 5 ) 


ind the equation for the progress variable 


c: 


dpc dpu{C _ dJk 


dt dxj dxk 


+ d> c 


( 6 ) 


where P, Tij and Jk are respectively the pressure, the viscous stress tensor and 
the molecular diffusive flux of c. Multiplying Eq. (5) by c and Eq. (6) by u,, then 



276 


D. Veynante & T. Poinsot 


adding and averaging the two resulting equations leads to a transport equation for 
pule. In a similar way, adding the averaged version of Eq. (5) multiplied by c and 
the averaged form of Eq. (6) multiplied by U{ leads to a transport equation for 
puiC. Subtracting the two resulting equations provides a transport equation for the 
turbulent flux pu'/c": 


(I) (II) (III) (IV) (V) 

+ pu”u> + c"F t v 


dUi 


(II) 


(HI) 


— &P 

„dp' 

"dJk 

r» dT ' k 

dxi 

dx i 

-“•0x t +C 

dx k 

(VI) 

(VII) 

(VIII) 

(IX) 


( 7 ) 


Bray et al.( 1981) studied each term in Eq. (7) and proposed some approxima- 
tions. For example, they explored the role of the mean pressure gradient term (VI) 
assuming that this term is so large that only cross dissipation terms (VIII and IX) 
can provide a balance, leading to a turbulent flux directly proportional to the pres- 
sure gradient. All terms in Eq. (7) may be extracted from our direct numerical 
simulations (Trouve et a/., 1994; Veynante et a/., 1995) to validate these analyses. 

Two main comments arise from Eq. (7). First, the mean pressure gradient appears 
explicitly in the source term (VI) ( c H dp/dxi ). Under the BML analysis, c" may be 
easily closed (Masuya & Libby, 1981): 


c n — c — c — 


1 + TC 


( 8 ) 


which is exact for an infinitely thin flame front. This quantity, being always positive, 
indicates a pressure decrease from fresh to burnt gases and tends to promote counter- 
gradient diffusion, or positive values of the turbulent flux pu”c. 

Another important feature concerning the present numerical simulations deals 
with volume and buoyancy forces. A constant volume force F” leads to a source 
term in Eq. (7) with a similar form for the pressure gradient term, whereas a 
constant acceleration force T, does not. However, the introduction of F? or has 
a direct influence on the mean pressure gradient dP/dx (term VI). So, terms (VI) 
and (XI) should be grouped to describe the effect of F? or I\. We start from the 
averaged momentum transport equation: 


dpui dpuiUj &pu' t 'u'j dP df ik 

dt ^ dxj dx 3 dxi + + 9 * + dx k 

For sufficiently large volume and/or buoyancy forces, an equilibrium between these 
forces and the mean pressure gradient can be used, which leads to the hydrostatic 
approximation: 

dP 
dx j 


F' + pTi 


( 10 ) 
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Therefore, to first order, term (VI + XI) becomes: 



Both Ff and T, induce a mean pressure gradient in the flow field (see Eq. 5). 
However, only I\ will directly affect the balance of pu"c" (Eq. 7). * Accordingly, 
our study of the influence of the mean pressure gradient on turbulent transport will 
be conducted using a constant acceleration I\. 

2.2 Physical and numerical issues related to pressure gradients in flames 

Theoretical models indicate that both normal and tangential pressure gradients 
influence turbulent flames. Masuya & Libby (1981) have studied confined oblique 
flames and have shown that, for a given pressure gradient, turbulent transport in 
normal and transverse directions are correlated. As a first step, we will only consider 
pressure gradients in the mean propagation direction x\\ only F? = F v and Ti = T 
may be nonzero. All pressure gradients are scaled by the pressure gradient inside 
the laminar flame zone given by: 

|VP, om | * p(s°,) 2 t/S° (12) 

where 6 {* is the unstrained laminar flame thickness that is obtained from the maxi- 
mum temperature gradient = Max( j^)/(Tt, — T u ). 

The pressure gradient VP/ am is created by dilatation inside the flame zone. It is 
large, but due to the thinnish flame front, the overall pressure jump between fresh 
and burnt gases remains small since A P/P ~ r'f(s° l /c) 2 where the ratio of flame to 
sound speeds is of the order of 0.001. On the contrary, volume forces or external 
pressure gradients are imposed on distances much larger than the flamelet thickness 
and will overcome the effect of dilatation in turbulent flame brushes. These effects 
and their relative importance in various flames may be quantified in terms of two 
quantities: the reduced external pressure gradient VP^ t , and the reduced mean 
pressure gradient VP+ eon . We consider these gradients positive when the pressure 
increases when going from the fresh to the burnt gases. Figures 1 and 2 illustrate 
how these gradients may be estimated in a turbulent flame brush of thickness /j 
with and without an external pressure gradient. With no external pressure gradient, 
corresponding to Fig. 1, the pressure jump across the flame will be conserved but 
spread over the flame brush thickness /(, so that the mean pressure gradient in the 
flame brush will be VP+ ea „ ~ VP,* m <5° /*&• In the case of an externally imposed 
pressure gradient, corresponding to Fig. 2, VP+ ean will be of the order of VP e * , 
everywhere in the flow; although the maximum instantaneous pressure gradient may 
still be found at the flame, VP e * t will be dominant over the flame brush thickness. 

Pressure gradients in real flames are imposed either by flame confinement, as 
in ducted flows, or by gravity. Typical values of VP,* m , VP* t , and VP* ean are 

* This result is due to the fact that, contrary to a constant acceleration, a constant volume force 
does not introduce buoyancy phenomena. 
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FIGURE 2. Pressure gradients in a flame with imposed pressure gradient. 

given in Table I for different turbulent premixed flames along with their dimensional 
values. In the case of “free” flames, the pressure gradient is imposed by the flame 
itself. The flame data correspond to the experiment of Schefer et al.( 1982). The 
flame brush thickness is of the order of the integral length scale l t = 1 cm. 


TABLE I. Typical pressure gradients for a propane/air flame (P=l atm, <f> = 1.2) 



s i 

m/s 

m 

T b /T u 

Pa/m 

|VP e * t | 

Pa/m 

IVP+,1 

|VP+ eon | 

Free flame 

0.32 

0.0005 

6.5 

250 

0 

0 . 

0.05 

1 g flame 

0.32 

0.0005 

6.5 

250 

~ 10 

0.5 

0.5 

Ducted flame 

0.32 

0.0005 

6.5 

250 

1000 

5 

5 











Effects of pressure gradients on premixed flames 


279 


These estimations indicate that the largest pressure gradients will be obtained 
in ducted flames where large pressure differences are created. Gravity alone will 
create smaller effects. At this point, it is worth discussing the differences between 
flames subjected to a constant acceleration, such as gravity, or to a constant pressure 
gradient. In the first case, the induced pressure gradient is pY , which is different 
in the fresh and burnt gases, see below. The induced pressure gradient will be 
constant in the latter case. For both cases, however, the pressure gradient will 
induce differential acceleration for fresh and burnt gas pockets in both cases thus 
leading to a modification of turbulent transport. Most authors therefore expect 
similar effects from constant acceleration or from pressure gradients. 

There are at least three ways to introduce pressure gradients in a direct numer- 
ical simulation of premixed turbulent flames: impose a constant volume force F v \ 
impose a pressure gradient through the boundary conditions; or impose a constant 
acceleration T, i.e. a volume force that is a function of the local density F g — pT, 
where T = g in the case of gravity. All of these techniques produce an imposed pres- 
sure gradient, see Appendix I. However, we have seen that the first solution leads 
to a flow where the pressure gradient is compensated everywhere by the volume 

force F v and has no effect on pu"c n . The second solution was investigated, but it 
is difficult to implement in a simulation if the mean flow remains one-dimensional, 
which is required for statistical purposes. In this paper, we will use only the third 
solution with various values of the acceleration T. 

3. Direct numerical simulation of premixed flames with pressure gradients 

The present direct numerical simulations were performed with a two-dimensional 
version of NTMIX. A complete description of this code may be found in Haworth 
&: Poinsot (1992) or Poinsot & Lele (1992). It solves the fully compressible Navier- 
Stokes equations with a single irreversible reaction Fuel — ► Products . Variable 
density as well as viscosity and transport coefficients are taken into account. The 
conservation equations solved by the simulation are: 


dp d , , 

at + a^ im) 


dpE 

~W 


dpui d 

sr + 

Q 

+ -fa\(P E + P) U i\ 


d{py) 

dt 




= 0 , 


j jpv I -p I 

= -W i +F - + fT ‘ + av 

+ i : (x W i ) + Q,i ' 


(13) 

(14) 


(15) 

(16) 


where 
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pE=\,,'E« 1 + p _ 

Z k=l 7 l 


(17) 

( dui duj 2 duk\ 


(18) 

Tij ~ ^ [dTj + dTi ~ 3 6ij dr k )’ 


w — wr/Yr = pYB exp (-13/ a) exp ^ 

-0(1-*) \ 
l-a{l-0)J‘ 

(19) 


In these expressions p is the mass density, p is the thermodynamic pressure, pE is 
the total energy density, Q designates the heat of reaction per unit mass of fresh 
mixture ( Q = —AhjYjt where Ah°j is the heat of reaction per unit mass of reactant). 
F t v and T , are the z- component of the constant volume force F v and the constant 
acceleration T, respectively. The reduced temperature is 0 = (T — T u )/(T b — T u ). 
The fresh gas temperature is T u , and T b is the adiabatic flame temperature for unity 
Lewis number. The activation temperature is T a . B is the pre-exponential factor 
and the coefficients a and /? are the temperature factor and the reduced activation 
energy, respectively: 

a = (T b - T u )/T 6 ; 0 = aT a /T b . (20) 

The mass fraction of the reactants Yr is nondimensionalized by the initial mass 
fraction of reactants Y# in the fresh gases: Y = Yr/Y£. This varies from 1 in the 
fresh gases to 0 in the burnt gases. 

We assume that the gas mixture is a perfect gas with constant molar mass and a 
specific heat ratio 7 = 1.4. The thermal conductivity A and the diffusion coefficient 
V are obtained from the viscosity coefficient p according to 

A = pCp/ P r and V = p/(pS c ), (21) 

where the Prandtl number P r and the Schmidt number Sc are constant. As a 
consequence the Lewis number L e = S C /P T is also constant. The viscosity p is a 
function of temperature: p = p u (T/T u ) b . 

The computational domain is L x by L y with N x by N y grid points. Two box 
sizes have been used: the small box has N x = N y = 257. The aspect ratio is L y IL x 
= 2.5. The second box, or large box, has N x — 257 and N y = 1025. The aspect 
ratio for this box is £ y /L r = 6.66. The following parameters have been used for 
both cases: 

TABLE II. Fixed parameters for direct numerical simulations of flames subjected 
to pressure gradients 


Re = cL x /v Le 

Pro 

Ta/T b 

T b /T u b 

s?/c 

6° /L x Nx 

12000. 1. 

0.75 

8. 

4. 0.76 

0.0159 

0.027 257 


where the speed of sound and kinematic viscosity in the fresh gases are denoted by 
c and v. 
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FIGURE 3. Pressure profiles in ID laminar flames without and with imposed 

acceleration g*. g* = 0 ( ); g* = —6.25 ( ); g* = 3.12 ( ) ; 

9 * = 6.25 ( ) 



Figure 4. Density profiles in ID laminar flames without and with imposed 

acceleration g*. g* = 0 ( ); g* = —6.25 ( ); g* = 3.12 ( ); g * = 6.25 

( ) 

4. Laminar flames submitted to pressure gradients 

First, one dimensioned lamineir flames eire computed without and with an imposed 
constant acceleration T. Introducing the reduced acceleration, g* = r<5 ; ° /(s® ) , 
which may be viewed as the inverse of a Froude number, four values of g* are 
considered here: g* = 0 (no imposed acceleration), g* = —6.25 (favorable pressure 
gradient), g* = 3.12 and g* = 6.25 (adverse pressure gradient). Pressure profiles 
are plotted as a function of the downstream locations in Fig. 3 for the four g* 
values. As expected, pressure gradients axe constant for each side of the flame front 
but decrease by a factor Tf,/T u = 4 between fresh and burnt gases due to density 
changes. The pressure drop due to thermal expansion is apparent for the g* = 0 
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FIGURE 5. Configuration for the numerical simulations, 
case. 

For all of the g* values considered, the laminar flame structure is not affected by 
the imposed acceleration: flame thickness, reaction rate and mass fraction profiles 
remain unchanged. Nevertheless, due to the pressure gradient, a weak change in 
density, similar to the one described in Appendix 1, is observed inside the fresh 
or burnt gases, see Fig. 4, but remains negligible compared to the one induced by 
thermal expansion. The same trend is noticed for the flow velocity which is modified 
by about 3% by the pressure gradient and a factor of 4 by thermal expansion. 

TABLE III. Numerical parameters for direct numerical simulations of 2D turbulent 

flames 

Case u'Js ? g* N x N y 


A 

5. 

3.5 

0 . 

257 

1025 

B 

5. 

3.5 

-6.25 

257 

1025 

C 

2. 

2.7 

0 . 

257 

257 

D 

2. 

2.7 

3.12 

257 

257 

E 

2. 

2.7 

6.25 

257 

257 


5. Turbulent flames submitted to pressure gradients 

The previously computed ID laminar flames are used as initial solutions for 2D 
flame-turbulence interaction simulations (Fig. 5). A Passot-Pouquet turbulence 
spectrum, with given turbulence intensity u' and integral length scale It, is super- 
posed on the combustion field (Haworth & Poinsot, 1992). Two sets of numerical 
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(a) 




(b) 


FIGURE 6. Superimposed instantaneous temperature and vorticity fields. Initial 
turbulence level u' 0 /s ? = 5. (a) No imposed pressure gradient ( g * = 0. - case A); 
(b) Favorable imposed pressure gradient (g* = —6.25 - case B). 


simulations have been conducted. The first set, consisting of runs A and B, starts 
from a high turbulence level ( Uq/$ *} = 5) and exhibits a classical gradient turbu- 
lence transport. The initial Bray number for this flow is Nb = 0.6. A favorable 
pressure gradient (i.e. VP+ ean < 0) is imposed with a positive acceleration ( T > 0) 
to reach a counter gradient diffusion situation. The second set, consisting of runs 
C-E, starts from a low turbulence level (uq/s® = 2), counter-gradient situation, 
with an initial Bray number of Nb = 1.5. Under an adverse pressure gradient (i.e. 
VPj^ean > 0), the flow is found to exhibit gradient turbulent transport. Numerical 
parameters are displayed on Table III where g* is the reduced imposed acceleration 
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FIGURE 7. Instantaneous pressure field. Initial turbulence level u' 0 /s = 5. (a) 
No imposed pressure gradient (g* = 0. - case A); (b) Favorable imposed pressure 
gradient (g* = —6.25 - case B). 


(g* = T8® /(s®) 2 ). The values of VP+ ean were chosen of the order of the pressure 
gradient found in the experiment of Shepherd et al.( 1982). 

5.1 Effect of the mean pressure gradient on the turbulent flame structure 

Instantaneous temperature and vorticity fields are displayed in Fig. 6 for an initial 
turbulence level u' Q /s ~ 5 without (case A) and with (case B) an imposed mean 
pressure gradient. Corresponding pressure fields are displayed in Fig. 7. The flame 
structures are quite different. Due to the favorable pressure gradient ( dP/dx <0), 
the wrinkling of the flame front is lower and the turbulent flame brush is thinner. 
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FIGURE 8. Superimposed instantaneous temperature and vorticity fields. Zoom 
from Fig. 6 Initial turbulence level u' 0 /s° = 5. (a) No imposed pressure gradient 
(g* = 0. - case A); (b) Favorable imposed pressure gradient ( g * — —6.25 - case B). 



FIGURE 9. Superimposed instantaneous temperature and vorticity fields. Initial 
turbulence level Uq/s° = 2. (a) No imposed pressure gradient (g* = 0. - case C); 
(b) Adverse imposed pressure gradient ( g * = —6.25 - case E). 
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Reduced time 

Figure 10. Reduced turbulent flame speed Sr/tf plotted as a function of the 

reduced time t/tf, where if — S^/s*} is a flame time, for different values of g *. 

(a) initial turbulent level Uq/sJ = 5 ; cases A ( ) and B ( ); (b) initial 

turbulent level u' 0 /s ° = 2 ; cases C ( ), D ( ) and E ( ). 

Despite similar minimum and maximum values, the pressure field is mainly domi- 
nated by vortices in case A whereas the pressure gradient, imposed by the constant 
acceleration T, is clearly apparent for case B. 

Close-ups of the temperature and vorticity fields of Fig. 6 are displayed in Fig. 8. 
As previously described, the flame front is less wrinkled in case C despite a simi- 
lar turbulence distribution in the fresh gases. The internal flame front structure, 
however, is strongly modified by the mean pressure gradient and buoyancy effects. 
Protrusions of high temperature levels are clearly apparent in fresh gases. Ac- 
cordingly, local flame front characteristics may be strongly affected, which will be 
investigated in the future. 

Instantaneous temperature and vorticity fields are displayed for cases C and E 
in Fig. 9. The initial turbulence level is lower (tto/sj = 2): without externally 
imposed pressure gradient, this flow exhibits a counter gradient turbulent transport 
as predicted by the Bray number criterion ( Nb = 1.5). In case E, an adverse 
pressure gradient is imposed and a transition towards gradient transport is expected. 
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FIGURE 11. Reduced turbulent flame thickness St/ 8° plotted as a function of 
the reduced time t/t f , where tf = 8°/s^ is a flame time, for different values of 

g* . (a) initial turbulent level Up/s® = 5 ; cases A ( ) and B ( ); (b) 

initial turbulent level u' 0 /s° = 2 ; cases C ( ), D ( ) and E ( ). 

The turbulent flame brush is determined from the maximum value of the c gradient 
(6 t = 1 /Max(dcfdx) 

The flame front wrinkling is somewhat increased by the adverse pressure gradient, 
due to the differential acceleration induced by buoyancy between fresh and burnt 
gases. 

5.2 Effect of the mean pressure gradient on global turbulent flame characteristics 

The global turbulent flame characteristics, namely the turbulent flame speed 
St and flame brush thickness St, are plotted in Figs. 10 and 11 as a function of 
reduced time for different values of g*. As expected^ from the previous flow-field 
visualizations, a favorable pressure gradient, i.e. dP/dx < 0, which is generally 
encountered in practical situations of confined turbulent flames, leads to a thinner 
turbulent flame brush and a lower turbulent flame speed. The decrease jn St may 
reach 30%. On the other hand, an adverse pressure gradient, i.e. dP/dx > 0, 
induces an increase in flame brush thickness and a higher turbulent flame speed. 
These results are in agreement with the influence of a constant acceleration T on 
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FIGURE 12. Transverse profiles of the turbulent flux pu"c" plotted as a function 
of the mean progress variable c for different reduced time step t + = t/tf where 
tf = <5® /s° is a flame time. Initial turbulence level Uq/s) 1 = 5. t + = 0.6 ( ); 0.9 

( ); 1.2 ( ); 1.5 ( ); 1.8 ( ). (a) No imposed pressure gradient 

(g* = 0. - case A); (b) favorable imposed pressure gradient (g* -■ —6.25 - case B). 

St theoretically predicted by Libby (1989). 

5.3 Turbulent transport pu"c" 

The transverse profiles of the turbulent flux pu"c" as a function of the mean 
progress variable c for the cases A and B at various times are shown in Fig. 12. 
Case A^without an imposed mean pressure gradient, is clearly of gradient type, 
i.e. ~pu"c" < 0, whereas the imposed favorable pressure gradient leads to a counter 
gradient turbulent transport. This finding is in agreement with the work of Bray et 
a/. (1982) and is expected from Eq. (7). Even in clearly counter gradient situations, 
the turbulent flux pu"c" is always negative, or of gradient type, at the leading edge 
of the flame brush, where c -» 0. As shown by Bray and his coworkers, these 
gradient zones allow flame stabilization. 

The total turbulent flux, i.e. pu'{c"dx, is plotted as a function of the reduced 
time for the different simulations in Fig. 13. Favorable pressure gradients promote 
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FIGURE 13. Total turbulent flux pu" c n dx plotted as a function of the reduced 

time t/tf where tf = 8® /s® is a flame time for the different cases, (a) initial 

turbulence level Uq/s*} = 5 : Cases A ( ) and B ( ); (b) initial turbulence 

level Uq/s ® —2 : Cases C ( ),D( ) and E ( ). 

counter-gradient diffusion, corresponding to a reduction of both the turbulent flame 
speed St and the turbulent flame thickness S T . On the other hand, adverse pressure 
gradients lead to an increase in St and St and induce gradient turbulent transport. 

5.4 Analysis of the pu n c n transport equation 

All terms in Eq. 7 may be obtained from direct numerical simulations. This anal- 
ysis serves to identify the dominant terms as well as the nature of their contribution. 
A typical direct numerical simulation evaluation of terms (I)-(X) appearing in the 
c-flux budget of Eq. 7 is presented in Fig. 14 for case C. The figure also displays the 
imbalance that was found when numerically closing the c-flux budget in Eq. 7. This 
imbalance is due to inherent numerical errors involved in the simulations as well as 
in the post-processing of the data. Its magnitude remains small, which suggests that 
the simulations can be used to analyze the variations of second-order moments. For 
instance, Fig. 14 shows that the dissipation terms (VIII) and (IX), which are gener- 
ally modeled together, are of the same order and act to promote gradient diffusion. 
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mean progress variable 


FIGURE 14. Variations of the different terms appearing in the c-flux budget across 
the turbulent flame brush. Case C: counter-gradient turbulent diffusion and no 
imposed pressure gradient. 

On the other hand, pressure terms (VI) and (VII), and the velocity-reaction rate 
correlation (X), strongly act to promote counter- gradient diffusion. The two source 
terms due to mean progress variable gradient (IV) and mean velocity gradient (V) 
tend to decrease the turbulent fluxes as expected and, accordingly, in the present 
counter-gradient situation act to promote gradient turbulent diffusion. The mean 
pressure gradient term corresponds to the pressure jump across the flame brush (see 
Fig. 17): 


dx ( ’dx 1 + re S T 


( 22 ) 


The fluctuating pressure term (VII) cannot be neglected as generally assumed in 
the models proposed to close the transport equation (7). 

A similar analysis is now performed for case E, where, due to the imposed adverse 
pressure gradient, the turbulent diffusion becomes gradient type, as indicated in 
Fig. 15. As expected, the mean pressure gradient term tends to promote gradient 
turbulent diffusion and corresponds to the imposed pressure gradient. Once again, 
the fluctuating pressure term (VII) is not negligible and acts to counterbalance 
the mean pressure gradient term (VI). In fact, the combined term (VI) + (VII) is 
mainly negative and corresponds to a gradient diffusion. The reaction term acts to 
promote counter-gradient diffusion. 

The budget of the transport equation for pu"c" for case B is presented in Fig. 16. 
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mean progress variable 

FIGURE 15. Variations of the different terms appearing in the c-flux budget 
across the turbulent flame brush. Case E: gradient turbulent diffusion induced by 
an adverse mean pressure gradient. 


The favorable mean pressure gradient acts to promote counter-gradient turbulent 
diffusion from term (VI). Once again, term (VII) tends to counterbalance term (VI). 
For such a situation, Bray et al( 1982) propose to neglect the pressure fluctuation 
effects (term VII) and assume that only the cross-dissipation terms (VIII 4* IX) 
can provide a balance to the mean pressure term (VI). In fact, from our simulation, 
the mean pressure term (VI) is balanced by the sum of the three contributions: the 
cross-dissipation term (VIII + IX), the pressure fluctuation term (VII) that cannot 
be neglected, and the source term due to gradients of c (IV). 

The mean pressure gradient across the flame brush may be simply modeled as 
the sum of two contributions: the imposed pressure gradient and the_pressure jump 
due to thermal heat release. As a result, the source term (VI) in pu n c n becomes: 



(c -% 


rss T 


g 1-c) 

1 + TC 


~ pr \ 


(23) 


where pT corresponds to the imposed pressure gradient. This expression is verified 
in the present simulations (Fig. 17). 


6. Theoretical analysis and modeling 

6.1 Model for turbulent flux without pressure gradient 

In this section, we first recall the derivation of a model for the turbulent flux of the 
mean progress variable c. Details about this modeling may be found in Veynante et 




292 


D. V eynante & T. Poinsot 



0.0 0.2 0.4 0.6 0.8 1.0 


mean progress variable 

FIGURE 16. Variations of the different terms appearing in the c-flux budget across 
the turbulent flame brush. Case B: counter-gradient turbulent diffusion induced by 
a favorable mean pressure gradient. 

a/.(1995). This derivation starts from a relation proposed by Bidaux & Bray (1994) 
connecting the flame front averaged fluctuating velocity, (u") a , to the conditioned 
unburnt and burnt mean velocities: 


«), = (c* - c)(«j 6 - u<«) 

where c* corresponds to the iso-c line used to define the flame location. 
From Eq. (4), the previous relation leads to: 


«>, = 


C( 1 - 2 ) *' 


This expression may be used to derive an estimate of u"c" via a model for the 
mean velocity fluctuation («"),; considering limiting cases of low turbulence lev- 
els where flow dynamic is mainly controlled by the thermal expansion across the 
flame brush, and high turbulence levels where the turbulent velocities dominate the 
flow induced by thermal expansion, Veynante et a/. (1995) proposed the following 
expression for (u" ) where index 1 corresponds to the direction normal to the flame: 


( u i)« = ( c * - <D( r5 l - 2m/) 


(26) 
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FIGURE 17. Comparison between direct numerical simulation data for cases B 

( ), C ( ) and E ( ) and modeled mean pressure gradient source 

term ( ) from Eq. (23) in the the c-flux transport equation (Eq. 7). 

leading to: 


w'/c" = c(l — c)(rsi — 2cm') (27) 

and to the Bray number criterion. Here u f denotes the rms velocity fluctuations and 
a is an efficiency function to take into account the low ability of small turbulent 
vortices to affect the flame front. 

6.2 Model for turbulent flux with pressure gradient 

Our objective in this section is to incorporate pressure gradient effects in the 
previous analysis. Pressure gradients induce differential buoyancy effects between 
cold heavy reactants and light hot products. We quantify this buoyancy effect 
through a characteristic velocity Ub (c), which is simply added to the two velocities 
used in Eq. (26): the velocity induced by the flame, (c* — c)rsi\ and the velocity 
induced by turbulence, 2(c* — c)au', 

( u i)» = ( c * “ c)(™L - 2au') + U B (c) (28) 

The estimation of Ub(c) is done as follows. In the fresh gases, J7 b(0) = U b B 
corresponds to the relative speed of a pocket of burnt gas (density diameter /). 
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Similarly a velocity J7 b(1) = Uq will be associated to the movement of unburnt gas 
pockets in the burnt products. Assuming an equilibrium between buoyancy and 
drag forces and using Stokes law for the drag of the pocket, U h Q is given by: 


U b B 


1 r / 2 
12 r + 1 v u 


r 


(29) 


where v u is the kinematic viscosity in the fresh gases. 

For a pocket of fresh gases with density p u and diameter /, the same analysis 
leads to a relative displacement speed [/£: 


Ul 


i— r 

12 v b 


(30) 


where v b is the kinematic viscosity in the burnt gases. 

Assuming a linear variation of the buoyancy velocity Ub(c) with c between the 
flame leading and trailing edges leads to: 


U B (c) 


1 r l 2 
12 r + 1 u u 


T [(l + (r + l) 1-n ) c — 


1 ] 


(31) 


where the kinematic viscosity ratio has been estimated as: 



and n = 1.76. 

For 3 < r < 6, we have 0.74 < 1/ (l + (r + l) 1_n ) < 0.8. In a first 
achieve simple expressions, the c-level retained to define the flame front 
chosen as c* = 1/ (l + (r + l) 1-n ), leading to: 



step, to 
may be 


U B (c 


1 T l 2 
12 c* (r + 1) v u 


F(c-c’) 


(32) 


Then, from Eq. (25) and (28), a simple model for the turbulent flux if'/ c" is: 


u'lc" - g(l -c)(rs L - 2W - $ — ( ;' + t) - r) (33) 

where a model constant is introduced to take into account the various limitations 
of the simplified analysis proposed here. 

The previous expression needs a characteristic length / corresponding to the typ- 
ical size of a pocket of fresh gases in a medium of burnt gases, and vice versa. As a 
first step, the integral length scale l t , a rough estimate of the flame front wrinkling 
scale, is used. Predictions from expression (33) axe compared with simulation data 
for the reduced time 1.8 in Fig. 18 using c* = 0.8 and /? = 0.12. The efficiency func- 
tion, a, is a function of the length scale ratio It/ Si and is obtained from previous 
direct numerical simulations (Veynante et a/., 1995) to be about 0.5 for the length 
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FIGURE 18. Comparison between direct numerical simulation data (Bold curves) 
and model predictions (thin curves) from Eq. 33 of the turbulent flux u"c" as a 

function of the mean progress variable c at reduced time 1.8 for cases A ( ), 

B ( ), C ( ), D ( ), and E ( ). a = 0.5, 0 = 0.12 and c* = 0.8. 

The characteristic length / is taken equal to the integral length scale l t 


scale ratio used here. The agreement between numerical data and model predic- 
tions is satisfactory. The influence of the imposed acceleration and the transition 
between gradient and counter-gradient transport are well predicted from Eq. (33). 

A simple criterion may be derived to predict the occurrence of counter-gradient 
turbulent diffusion, i.e. u"c" > 0), from Eq. 33: 




0 


/ 2 r 

12c* (r + 1)u u sl 


> 1 


(34) 


The effect of the constant acceleration is to introduce a coefficient for the Bray 
number Nb defined by Veynante et a/. (1995), see Eq. (3). Ng may be rewritten as: 




1-0 


12 c* (r + 1) 



> 1 


(35) 


where Ref = Sisi/u u is a flame Reynolds number. 
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This modified Bray number may be estimated from our numerical simulation 
using the initial turbulence values (/ = l tj v! — Uq, sl — s 1 /), a = 0.5 (from 
(Veynante et a/., 1995), f3 = 0.12 and c* = 0.8. Results axe summarized and 
validated in Table IV; it appears that > 1 flows indeed exhibit counter-gradient 
diffusion. 

The same analysis may be easily extended to the case of an externally imposed 
pressure gradient dP /dx, leading to 


ub ® * (£) 1(1 + <t + i), " n) i] 


u"c" = c(l — c)( tsl — 2 au' — 


1 “'V \2p u v u c* \dx ) J 

This expression leads to the following criterion for counter-gradient turbulent trans- 
port: 

= (38 > 



Si 

PuTsf 


which is the reduced pressure gradient (i.e. the pressure gradient made nondimen- 
sioncdized by the pressure gradient across the corresponding laminar flame). The 
previous result is based on a simple analysis, assuming an equilibrium between 
buoyancy and drag forces for a pocket of fresh (burnt) gases in a medium of burnt 
(fresh) gases. This approach is probably too crude and has to be improved to de- 
scribe the instability mechanis ms o f an interface between two gases with different 
densities. The dependence of u n c n with the square of the integral length scale l t 
needs to be examined and validated, especially for large values of l t . 


TABLE IV. Estimation of Ng. GD (CGD) refers to gradient (counter-gradient) 
turbulent diffusion 


Cast 

t w 

UlSi 

9* 

GD/CGD 

N T b 

A 

5. 

3.5 

0 . 

GD 

0.60 

B 

5. 

3.5 

-6.25 

CGD 

1.3 

C 

2. 

2.7 

0 . 

CGD 

1.5 

D 

2. 

2.7 

3.12 

GD 

0.96 

E 

2. 

2.7 

6.25 

GD 

0.43 


In the experiment of Shepherd et al.( 1982), u' « 1 m/s and /t ~ 1 cm, leading 
to Bray numbers Nq « 0.9 and N V B « 20. Due to the large length scale ratio, 
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l t /6i « 20, a = 1 , (cf. Veynante et aZ.1995). These results axe compared in Table 
V with the experimental estimation of the turbulent flux pu'/c". Accordingly, the 
unconfined case without pressure gradients lies in the transition region between 
gradient and counter-gradient diffusion. In fact, counter-gradient diffusion is found 
in the experiment, but the turbulent flux remains quite low. On the other hand, 
the confined case with an externally imposed pressure gradient, N B measures the 

intensity and the sign of the turbulent flux pu”c". The large value of N V B for the 
confined experiment is questionable because of the large length scale ratio ( It /Si ). 

TABLE V. Estimation of N B in the Shepherd et ai.(1982) experiment for the con- 
fined (with pressure gradients) and unconfined (without pressure pressure gradi- 
ents). Experimental estimations of pu"c"/pUo , where Uq is the reference burner 
inlet velocity, are provided, u 1 « 1 m/s, l t « 1 cm. 


Case 

K 

pu"c"/pUo 

unconfined 

0.9 

0.0043 

confined 

20. 

0.051 


7. Conclusion 

The influence of a constant acceleration T on a turbulent premixed flame is stud- 
ied by direct numerical simulation. This acceleration T induces a mean pressure 
gradient across the flame brush, leading to a modification of the turbulent flame 
structure due to differential buoyancy mechanisms between heavy cold fresh and 
light hot burnt gases. Such a pressure gradient may be encountered in practical 
applications in ducted flames. 

A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt 
gases, is found to decrease the flame wrinkling, the flame brush thickness, and the 
turbulent flame speed. A favorable pressure gradient also promotes counter-gradient 
turbulent transport. On the other hand, adverse pressure gradients tend to increase 
the flame brush thickness and turbulent flame speed, and promote classical gradient 
turbulent transport. 

The balance equation for the turbulent flux pu n c n of the Favre averaged progress 
variable c is also analyzed. The first results show that the fluctuating pressure 
term, (c n dp f dx), cannot be neglected as generally assumed in models. Simple mod- 
els assuming that a high mean pressure gradient may only be balanced by the 
cross-dissipation term seem too approximate. This analysis has to be continued 
to compare simulation data and closure schemes proposed for the ~pu"c n transport 
equation. 

The analysis developed by Veynante et a/.(1995) has been extended to imposed 
acceleration and mean pressure gradients. A simple model for the turbulent flux 
u"c n is proposed and validated from simulation data. Then, a modified criterion is 
derived to delineate between counter-gradient and gradient turbulent diffusion. In 
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fact, counter- gradient diffusion may occur in most practical applications, especially 
for ducted flames. 
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Appendix I: tests of NTMIX with constant force or accelerations 

To validate the implementation of volume forces in the simulation, simple exam- 
ples of one-dimensional flows with constant force or acceleration have been run. We 
consider an isentropic one-dimensional flow submitted to a force F v . F v is either 
constant or equal to pj. The governing equations are: 

djpu) _ Q 

dx 

d(pu 2 + p ) _ F 

dx 

For constant volume force F vy F is constant and equal to F v and for constant 
acceleration, F = pT. 

This system is integrated once to give: 


pu = po^o 

du «o F 

di = "^> 0 /u) 7+1 (i. -M*(u/u 0 y+') 

where index 0 designates the inlet condition. M = u/c is the local Mach number 
and 7 = 1.4. 

If the local Mach number M is small, this system may be integrated easily to 
give: 

For constant volume force F v : 

u(x) = u 0 (i - y lh 
■*0 

p(x) = p 

x'o 

P(x) = Po - F V X 

Note that = —F v 
For constant acceleration T: 

u(x) = « 0 (1 + 

c o 

eW = «(i + i2± ^) ,/h+11 

c o 
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p( X ) = p 0 ( i + (7+ 2 1)ri r /(7+1) 

c o 

where Co is the inlet sound speed (co = jPo/po)- 
For small values of T, = poT 

Fig. 19 shows (for a constant volume force F v /pqCq = 0.01) that the simulation 
results match these analytical expressions well. (A similar agreement is obtained in 
the case of a constant acceleration) 



FIGURE 19. Tests of code for one- dimensional flow with constant volume force 
Fv/poCo = 0.01. Comparison between simulation data: velocity (■ ), pressure (o ), 
density (* ) and analytical results ( ). 



